Modification of spin mixing of spinor BEC by cavity QED coupling 
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Dressed states of spinor Bose- Einstein condensates of spin-1 atoms coupling with optical cavity 
modes with far off resonance frequency are investigated. The exact solution of time evolution 
of population of spin component is derived, and the numerical result shows that the evolution is 
different from spin mixing without the coupling. Due to the coupling with the atoms, the photon 
state also evolute to different optical cavity modes. 
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I. INTRODUCTION 

Cavity Quantum Electrodynamics (cavity QED) [l], Q 
describes the coupling between atom and confined cavity 
mode of electromagnetic field by coherent interaction. It 
gives a great chance to reveal the nature of coherent in- 
teraction between matter and field, and to various appli- 
cation such as quantum information and quantum state 
engineering. The atom can be described by two-level sys- 
tem, and the cavity mode is described as quantized light 
field. The coupling between two-level system and quan- 
tized light field is described by the Jaynes-Cummings 
model [3(, which describe the coupling strength as the 
speed of energy transfer between two subsystem (here 
the excitation energy of the atom is close to energy level 
of cavity mode, which is called near resonance). When 
the energy level distribution of the atom is more compli- 
cated, three-level system or other multi-level system is 
needed to describe the atom. But the description of cou- 
pling is similar. Some multi-level system give interesting 
phenomenon such as electromagnetical induced transpar- 
ent Q. 

When there are more than one atom,the system can be 
described by Tavis-Cummings model _5|, where N atoms 
are assumed to couple to the same cavity mode with iden- 
tical coupling strength. This lead to enhance of the cou- 
pling strength to >/~N time. When the temperature of the 
system is lower than critical temperature, the N atoms 
couple with each other coherently and form Bosc-Einstein 
condensate (BEC). In this case, all atoms have the same 
quantum state, which ensure the same coupling strength 
with cavity mode. As a result, the Tavis-Cummings 
model is very appropriate for this system. A simpler 
model is to describe the BEC of the N atoms by mean 
field theory. A dressed BEC state @ and a dark state 
Q is investigated by the mean field theory description. 
When the energy level of cavity mode is much different 
from excitation energy of the atoms (far off resonance), 
the coupling can be described by effective Hamiltonian 
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Q , which is given by solving the Schrodinger equation of 
single atom and cavity mode. The dynamic property of 
dress BEC state of far off resonance is also investigated 
under mean field theory method [f|. 

When the N atoms is trapped in an total optical trap, 
the spin freedom of the atoms is librated. When temper- 
ature is low enough, spinor BEC is form which has 
much potential in application of quantum computer and 
quantum information [lfj | . Spinor BEC can be described 
by either mean field theory or analytical express of the 
quantum state. A former work about fractional parent- 
age coefficients (FPCs) gives a chance to describe spinor 
BEC of spin-1 atoms in a simpler way [TJ O ne of 
the most interesting phenomenon of spinor BEC is spin 
mixing, where the average atomic number of each spin 
component evolute versus time [l3l [Tij . The pattern of 
time evolution is found to be decided by symmetric prop- 
erty and initial condition of the system [la ] . If the spinor 
BEC is placed inside a far off resonant cavity, the atoms 
will couple to the cavity mode, which change the symmet- 
ric of the system as well as the effective coupling strength 
between atoms. As a result, the evolution of spin com- 
ponent is supposed to be modified. In this paper, we will 
deduce the effective Hamiltonian to describe the far off 
resonant coupling of spinor BEC and cavity mode, and 
figure the dressed state of spinor BEC which is eigen- 
state of the Hamiltonian. And then we investigate the 
modification of spin mixing by the coupling. The spinor 
BEC system is described by FPCs, which can give an 
analytical solution to the problem [1J 



II. DESCRIPTION OF DRESSED SPINOR BEC 
STATE 



The system we are going to investigate consist of ./V 
87 Rb atoms which are optically trapped to form a spinor 
BEC, and located inside a high-Q optical cavity. As- 
sume that there are only two optical modes of frequency 
lu Ci which are left and right circular polarization opti- 
cal mode in the cavity. The two optical modes have 
opposite angular momentum H along z direction, which 
are indexed by fi = ±1. The ground state of 87 Rb is 
5 2 Si/ 2 ,F — 1, and the lowest energy state of excited 
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state of Dl line is 5 2 P 1/2 ,F = 1 [13. When lo c is 
smaller than and at the same scale of the frequency of 
5 2 Si/ 2 ,F = 1 — > 5 2 Pi/ 2 ,-F = 1 transition w e , we can 
neglect other transition. As a result, the non-coupling 
Hamiltonian of the N atoms with spin dependent inter- 
action with the optical frequency removed is 

N 2 

H ° = E(|^ + C/ ^)+ E fti\ei,<r)(ei,a\) 

i=l <r=l,0,-l 

+ E^ w 

where is momentum of the ith atom, m is atomic mass, 
U(ri) is potential of external optical trap; S = cu e — u> c , 
\ei,a) is the 5 2 Pi/2, F = 1 state of the ith atom with spin 
component a (\gi, a) is the 5 2 <Si/2, F = 1 state of the ith 
atom with spin component a); Vi : j = <5(rj — rj)(co+C2Fj- 
Fj) is spin dependent interaction between atoms, with F 
being spin operator and cq 2 being interaction constant 

a-. 

Since the Q factor of the optical cavity is very high, and 
temperature is very close to zero, the decay of photon and 
non-coherent interaction between atoms and photon can 
be neglected. According to the selection rule, if the atom 
absorb (emit) one photon from(to) the cavity, the allowed 
transitions have Act = ±1. Therefore, the Hamiltonian 
of coupling between atom and cavity mode is given by 

N 

H c = X! M V|e 4 ,0)(. 94 ,-l|ai + \e l ,l)(g l ,0\a 1 

i=l 

+|e i ,-l)<0 i ,O|a_i + \ei,0)(gi,l\a-i + c.c) (2) 

where and ci+ are annihilation and creation operators 
of optical modes of polarization /1 = ±1, CIq = dE/H is 
the strength of the atom-field coupling, d is the atomic 
dipole-matrix element, and E — y/fvjj c /2eoV is electric 
field of single photon state, V being mode volume. In far 
off resonance region, we have 5 >> fin. 

Using similar process in reference [6(, we can get the 
effective Hamiltonian of the coupling system without in- 
volving the atomic excited state 

N 2 

^ = Et^ + ^o 

i=l 
Ml 2 

+ —g-(-\9i,-l){9i,-l\aiai ~ \9i, l|ai x a_i 
+ \gi,-l)(giA\afa-i + \giA)(9i,-l\a-iai)} 

+ E^- ( 3 ) 

It is obvious that the last two single particle terms are 
given by the process that the atom absorb one photon 
from one mode and emit one photon to the another mode. 
Note that the Hamiltonian ([3]) conserve the total angu- 
lar momentum at z direction M.H and total number of 



photons n = ni + n_i. In order to obtain the eigen state 
of the Hamiltonian, we first give the eigen state of the 
Hamiltonian without coupling between atom and optical 
mode (f2o — 0), which is direct product of eigen state of 
the N atoms and eigen state of optical mode. And then 
calculate the matrix element of the Hamiltonian with 
coupling (Oq 7^ 0) under the basic states given by the 
direct product states. The eigenstate of the Hamiltonian 
is obtained by diagonalizing the matrix. 
The direct product states is given as 

|5,M;ni,n_i) = \S, M) N ® |m,n_i) (4) 

where \S,M)u is eigenstate of the N atoms with total 
spin S and z component of total spin M, |ni,n_i) is 
optical state with n\ photons in /i = 1(— 1) po- 

larization state. The z direction angular momentum of 
the whole system is M.H = (M + m — n-i)%, and eigen 
energy of the states is Es, ni ,n-i = S(S + 1)J, where J 
is a constant decided by spin dependent interaction C2. 
The atomic eigenstate is given as 

N 

\S,M) N = #l M l[ ( t>(r i ) (5) 

i=l 

where is the total symmetric spin state of N atoms, 
</>(r) is the single particle spatial wave function. Note 
that we use second quantized Hamiltonian to describe 
optical mode, while use first quantized Hamiltonian to 
describe the N atoms. The total symmetric property of 
the N Bosonic atoms is included in the total symmetric 
atomic wave function. The total symmetric spin state 
can be expanded as 

*s, M = «5 [x(^s+i\s,m + b% (6) 

where a% = {[1 + (-1) N ~ S ](N - S)(S + 1)/[2N(2S + 
1)]}V2 and b% = {[l + (-l) N - s ]S(N + S+l)/[2N(2S + 
l)]} 1 / 2 are FPCs given in reference [Hj], x{i) is s P m state 
of the ith spin-1 ground state atom, which is couple to 
total symmetric spin state of the other TV — 1 atoms by 
Clebsch-Gordan coefficient. 

The eigenstates with quantum number n) can 

be expanded by the direct product states ^ with the 
same value of Mh and n = n\ + n_x- Since we are 
interested in the quantum effect of the system, we focus 
on one photon case, i.e. ri\ + n_i = 1, but the extension 
to more photon case is straight forward. When there 
is only one photon, only the direct product states with 
quantum number [M = M. — = l,n_i = 0) or 

(M = M, + l,ni =0, n_i = 1) can be basic states. As a 
result, the eigenstate can be expressed as 

\M,v)= cs,m-i\S,M~1;1,0) 

S>=M-l 

+ E c» SM+1 \S,M + 1;0,1) (7) 

S>=A4 + 1 
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where v is the index of eigenstate, c v s M±1 is expansion 
coefficients, and the summation satisfy N — S is even. 
Making use of equation ([6|), the matrix elements are ob- 
tained as 



(S,M - m; {fi}\H eff \S,M - m; M) = S(S + 1) J 

-^Ul M ^f + {Bl M _^f] (8) 



(S + 2 ) M- f i;{(i}\H ef f\S,M- f x;{n}) = 



m^n aN 



,B 



N 



S,M-ix,-ix JJ S+2,M—ti,—n 



(9) 



<5-2,X- M ;{ / i}| J ff e //|5,X- Ai ;{/i}) = 

hQ o N A N r" rim 



hQ 2 N 

(S, m + m; {-M}|# e //|S, - m; M) = -f- 



(^-S,M+tJ-,^S,M-^,-^ + ^S,M+fJ-,^S,M-iJ,-ij) (11) 



iJV 



A' 



(5 + 2, X + fi; {-iM}\H eff \S, M - fi; {»}) 



HQqN n 



(S - 2, M + ,i; {-n}\H eff \S, M~m {n}) 



m^n aN 



B 



N 



S-2,M+ij,,ij, "S.M- fj,,-fj, 



(13) 



where {fi = 1(— 1)} stand for optical state |1,0) 
( 1 0, 1> ) , -A-s,M,n = a S^S+l,l,M-n,ij, an d Bg M/i = 

h s C s-i,i,M-w with C k%MuM 2 bein g Clebsch- 
Gordan coefficient The eigenstate \M.,v) is the 

dressed spinor BEC state. 



III. MODIFIED SPIN MIXING 

In experiment of spin mixing, the initial atomic state 
is a Fock state |iVo, M) that each spin component has fix 
number of atoms. Nq is atomic number of spin compo- 
nent 0, Mh is z direction angular momentum. Atomic 
number of spin component 1 and -1 can be expressed 
by No and M. There are two processes that make the 
system evolute to the other Fock states. Process one is 
scattering between a pair of atoms, which can makes the 
spin component 1 and -1 jump to and 0, and vice versa, 
due to spin dependent interaction. The system evolute to 
Fock state with the same M. Process two is absorbtion 
of a /i = 1(— 1) photon and then emission of a fi = —1(1) 
photon by an atom. The system evolute to Fock state 
with M ± 2, i.e. the z direction angular momentum of 
the spinor BEC increase (decrease) 2H. 



Assume that the initial photon state is |1,0), i.e. the 
initial state of the system is \Nq, M; 1,0). The time evo- 
lution of wave function is given by 

|¥(f)) = e- iH -f^ n \N o ,M;l,0) 

= E E E \ S '' M - w {ti)@',M - n {n}\e-* H 'f^ h 

S',n v S 

\M, u)(M,v\S, M; {fi })(S, M\ M\N , M; {mo}) {mo =i } 
= ]T \S',M m {M}><^, {M0=1} (*) (14) 

S',M 

where Ai = M + 1. The matrix element between Fock 
state and total spin eigenstate (S, M — 1\N ,M — 1) is 
given in reference [HI]. After extracting the spin state 
of the Nth atom from the total spin eigenstate \S',A4 — 
H; {/x}) with equation ((6]), we can calculate the evolution 
of population of each spin component as 



p - pn=i _i_ DM=-1 



(15) 



S 



(16) 



where N — S is even and N—S'\s odd. The probability of 
finding the photon in cavity mode {/z} is given as P£ Q = 
S CT =±i o a' 11 l ne initial photon state is |0, 1), the 
same equation with M. = M—l and {/io = —1} describes 
the time evolution. If the initial photon state stay in two 
cavity modes with the same probability, the evolution is 
given by P a = (P ft0= i, CT + P M0= _i )O -)/2. We can see that 
the time evolution is decided by symmetric of the atomic 
system as well as coupling strength and spin dependent 
interaction strength. 

Numerical result of 300 Rb s7 atoms with various initial 
condition is shown in figure [TJ J is estimated by Thomas- 
Fermi approximation as J = h(ujf rap /N) 3 / 5 /X with 
X = 1.52 x 10 4 and w trQp = 514.5-ffz being frequency 
of the optical trap [13J, iv e = 2vr x 377.107TPz = 2w c , 
and Qq = 2tt x 1MHz which is available in experiment 
2]. The time is normalized by t per i OC i — nH/J, which 
is period of time evolution when there is no coupling be- 
tween atoms and cavity. From figure [1] (a), comparison of 
(i) and (ii) shows that time evolution of coupling system 
has oscillation at time zone that non-coupling system is 
stationary. In figure (iii), the time evolution of proba- 
bility of photon state at each polarization mode becomes 
random and fast oscillating around 0.5 after a short time 
of smooth relation. Changing to another initial condition 
in figure (b), time evolution of coupling and non-coupling 
system are almost the same. The photon state keeps os- 
cillation between two polarization modes and relax to be 
stable. In figure (c), it was predicted in a former work 
that this kind of initial condition result in a long stable 
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evolution and then a short plus at time t per i d/&. Figure 

(ii) shows that the coupling makes the stable evolution 
contains a small amplitude oscillation. The frequency of 
the oscillation small at early time and keep increasing, 
and then become a multi-plus at time t per i 0( i/8. In figure 

(iii) , it is shown that the photon state oscillates at a very 
large frequency between two polarization modes with a 
quasi-constant amplitude. 

IV. CONCLUSION AND DISCUSSION 

We have model the coupling system of TV Rb sr atoms 
and high Q optical cavity mode by two level approx- 
imation, i.e. only the transition 5 2 Si/ 2 ,F = 1 — > 
5 2 Pi/2,F — 1 of D line is consider. A more accurate 
model could be including other transition in D line of 
Rb 87 . Because the frequency of optical cavity mode is 
far off resonance from transition of D line, the opera- 
tor of excited state in the Hamiltonian can be removed, 
which result in an effective Hamiltonian. The effective 
Hamiltonian describe the coupling by two photons pro- 
cess, absorbtion of one photon from one cavity mode 
and then emission of one photon to another cavity mode. 
The eigenstate, or the dressed spinor BEC state, is ob- 
tained by diagonalizing the Hamiltonian under the basic 
set of direct product states, which are consist of total 



spin eigenstate and photonic Fock state of non-coupling 
Hamiltonian. The FPCs exhibit strong power in deduc- 
ing the matrix elements of the Hamiltonian. With the 
dressed states in hand, the time evolution of wave func- 
tion and average population in each spin component is 
derived. From the evolution equation (fT4")) (fT5|) and (|16p , 
when the system evolute to the direct product state with 
{/i =/= Mo}, the photon evolute from the initial cavity mode 
to another cavity mode, and at the same time the spinor 
BEC evolute to the state with different z direction an- 
gular momentum. This process result in entanglement of 
the cavity modes and the spinor BEC. The numerical re- 
sult show that the time evolution with some special initial 
condition is greatly modified by the cavity QED coupling. 
This effect is supposed to be observable experimentally. 
The time evolution with some other special initial con- 
dition, the oscillation of probability of finding a photon 
in each cavity mode is regular, which could be useful in 
quantum entanglement application or exploration. 
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i M=0, N Q =N-M (b) M=N/2, N Q =N-M (c) M=0, N Q =(N-M)/2 
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FIG. 1: Column (a) (b) and (c) is numerical result with differ- 
ent initial condition marked at the subtitle. Atomic number 
is N = 300. Row (i) and (ii) is time evolution of the aver- 
age population at spin component a = O(black), l(red), and 
-l(blue) without and with coupling between atoms and cavity, 
respectively. Row (iii) is time evolution of probability of the 
photon state at fi = l(blue) and -l(red) polarization mode. 



